Exact Quantum State of Collapse and Black Hole Radiation 
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I. INTRODUCTION 

The end state of the gravitational collapse of a compact 
object is an important problem for which a full under- 
standing has not yet been achieved. It is not completely 
clear what changes a quantum theory of gravity will bring 
about in our understanding of Hawking radiation, the en- 
tropy of a black hole, the information loss paradox and 
the nature of the gravitational singularity. This is be- 
cause it has not been easy to develop dynamical models 
of quantum gravitational collapse in the various candi- 
date theories of quantum gravity. 

Over the last few years, useful and interesting devel- 
opments in formulating the quantum gravitational col- 
lapse problem have taken place, within the framework 
of canonical quantum general relativity. Kastrup and 
Thiemann Q and Kuchaf 0, E| have pioneered the 
technique of midisuperspace quantization, a canonical 
quantization of a gravitational system with symmetries, 
and yet possessing infinitely many degrees of freedom. A 
prime example of this is the canonical quantization of the 
Schwarzschild geometry associated with an eternal black 
hole 0. 

We have applied these earlier works to develop a canon- 
ical description [f| of the collapse of a spherical time-like 
dust cloud, which is described by the LeMaitre-Tolman- 
Bondi family of metrics |(J. Quantization of the classi- 
cal dust system leads to the Wheeler-DeWitt equation 
for the wave-functional describing the quantum collapse. 
The quantum Schwarzschild black hole then becomes a 
special case of our model. We have successfully employed 
this equation to derive the Bekenstein mass spectrum 
and statistical entropy of the charged and uncharged 
black hole [8| and, in a recent work we showed that 
the semiclassical (WKB) treatment of the Schwarzschild 
black hole in this canonical picture describes Hawking ra- 
diation. The midisuperspace program has also been used 
to describe the collapse of null dust, with interesting con- 
clusions fijl \12L Ha. Ilj, [jjl . All these works demonstrate 



the overall consistency of the program. 

Apart from the usual operator ordering ambiguities, a 
problem with the midisuperspace quantization program 
is that the Wheeler-DeWitt equation involves two func- 
tional derivatives taken at the same point and requires 
regularization in all but the effectively minisuperspace 
models. Any attempt at regularization, however, leads 
to ambiguities that are not easy to resolve. Our pur- 
pose, in the present article, is to address this problem by 
applying a lattice regularization scheme to the canonical 
quantization of LeMaitre-Tolman-Bondi (LTB) models. 

We first construct a stationary state solution of the dif- 
feomorphism constraint. Upon discretization, the wave- 
functional reduces to a tensor product over shell wave- 
functions, one at each lattice point. The Wheeler-DeWitt 
equation is well defined on the lattice as an infinite set of 
Schroedinger equations, one for each shell. The results 
are independent of how the lattice is discretized when the 
size of all lattice cells approaches zero. Even though the 
regularization is on a lattice, diffeomorphism invariance is 
recovered when the lattice spacing approaches zero. The 
shell wave-functions reproduce Hawking's thermal radi- 
ation near the horizon in the WKB limit and the first 
quantum gravitational correction to Hawking's formula 
renders the spectrum non-thermal. We believe that our 
construction opens up a useful avenue for investigating 
some of the important unresolved issues related to black 
hole physics mentioned above. 

The paper is organized as follows: in section [H] we re- 
view our pervious quantization of the LeMaitre-Tolman- 
Bondi collapse models. In section IIIII we construct a 
lattice to regularize the Wheeler-DeWitt equation and 
discuss its properties. Solutions of the Wheeler-DeWitt 
equation are discussed in section Hvl and Hawking radia- 
tion is reexamined in section We conclude with a few 
comments in section IVll 
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II. CANONICAL QUANTIZATION OF LTB 
COLLAPSE 



The epoch R = describes a physical singularity, whose 
singularity curve 



The LeMaitre-Tolman-Bondi models .6] constitute a 
complete solution of the Einstein equations for a matter 
continuum of inhomogeneous dust, i.e., they are solutions 
of the spherically symmetric Einstein's field equations, 
G^ v = —SirGT^v, with vanishing cosmological constant 
and with stress-energy describing inhomogeneous, pres- 
sureless dust given by T^ v = eu^u v . The solutions are 
determined by two arbitrary functions, F(r) and E(r), 
of a spatial label coordinate, r, and given in co-moving 
coordinates as 



ds 2 = dT 2 

F 



R'- 



l + E 



-dr 2 - 



R 2 R 



u E ' Tr 



(i) 



(we have set 8irG = 1 = c) where R(t, r) is the physical 
radius of a spherical shell labeled by r. A prime rep- 
resents a derivative w.r.t. r and a star (*) represents a 
derivative with respect to the dust proper time, r. The 
arbitrary functions, E(r) and F(r), are interpreted re- 
spectively as the energy and mass functions. The en- 
ergy density of the collapsing matter is e(r, r), and the 
negative sign in the third equation above is required to 
describe a collapsing cloud. Its general solution is given 
up to an arbitrary function ip( r ) 01 the shell label coor- 
dinate. This arbitrariness reflects only a freedom in our 
choice of units i.e., at any given time, say t , the function 
R(t , r) can be chosen to be an arbitrary function of r. 

The mass function, F(r), represents the weighted mass 
(weighted by the factor yl + E) contained within the 
matter shell labeled by r. If a scaling is chosen so that the 
physical radius coincides with the shell label coordinate, 
r, at r = 0, then it can be expressed in terms of the 
energy density at r = according to 



F(r) = J e(0,r)r 



l dr, 



(2) 



while the energy function, E(r), can be expressed in 
terms of the initial velocity profile, v(r) — R*(0,r), ac- 
cording to 



E(r) =w 2 (r)- - / <(0.r)r<ir. 



(3) 



showing that it represents the total energy (gravita- 
tional plus potential) of the matter-gravity system. The 
marginally bound models, which we will consider in this 
paper, are defined by E(r) — 0. For the scaling referred 
to above, we must choose ip( r ) = r5 > whence the solution 
of can be written as 



RHr,r) 



(4) 



r(r) 



2ri 



(5) 



gives the proper time when successive shells meet the 
central physical singularity. Various models are obtained 
from choices of the mass function, F(r). For example, the 
Schwarzschild black hole is the marginally bound solution 
with F(r) = 2M, a constant. 

Reinserting the constants c, h and G, and defining the 
dimensionless variable x — r/l p (where l p = y^hG/c 3 
is the Planck length) and the dimensionless mass func- 
tion f{x) = F/l p , the wave-functional of the quantized 
Einstein-Dust system can be shown p| to satisfy the 
Wheeler-DeWitt equation, 
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r 



5t 2 {x) 5Rl[x) Al 2 \T 



*[T,R,f]=0, (6) 



where the prime denotes a derivative with respect to x. 
The upper sign refers to the region outside the horizon, 
R > F, and the lower sign to the region inside, R < F. 
We have also used T = 1 — F/R and a radius function 



R* 




e (—oo, +oo), 



(7) 



defined so that the DeWitt configuration space metric is 
manifestly flat. Invariance under spatial diffeomorphisms 
is implemented by the momentum constraint 0, 



RL 



SR, 



Sf 



*[t,R, /]=0. 



(8) 



To complete the quantum theory, one must define an in- 
ner product on the Hilbert space of wave-functionals. In 
, we defined it in a natural way, by exploiting the fact 
that the DeWitt super-metric yi| is manifestly flat in the 
configuration space (r, i?*), in terms of the functional in- 
tegral 



fl,(0) 



(9) 



Equations JJjJ EJ clearly imply a specific choice, albeit a 
natural one, of operator ordering, for which no justifica- 
tion deeper than the fact that they were shown to repro- 
duce the Bekenstein mass spectrum 0, of the eternal 
black hole and its Hawking radiation in the WKB ap- 
proximation [3] can be given at this time. This defect 
is unfortunately endemic to current theories of quantum 
gravity. 

Any functional that is a spatial scalar will obey the 
momentum constraint. In particular the functional 
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%{t,RJ] = exp 



ic 
21 



dxf'(x)W(r(x),R(x),f(x)) 



p Jo 



(10) 



is a spatial scalar if W(t, R, F) has no explicit depen- 
dence on x (equivalently, r) and we will use this as a 
solution ansatz. If we further require i|l(J|) to represent a 
stationary state then, noting that F'(r)/2 is the proper 
energy density of the collapsing cloud, we choose 



III. 



LATTICE CONSTRUCTION 



W(t,RJ) = t+U(RJ). 



(11) 



In a classical collapse problem, f(x) is generally mono- 
tonically increasing up to the boundary of the star, which 
is specified by some label Xb- The exterior of the star is 
then matched to the Schwarzschild exterior and f(x) is 
taken to be constant for x > x^. While this is a reason- 
able condition to impose upon the mass function on the 
classical level, it is too strong in the quantum theory. To 
account for Hawking radiation, it is necessary for f(x) 
to be a monotonically increasing function of x over the 
entire range of x, as we do not expect to have a sharp 
boundary between the collapsing matter and its exterior 
owing to the evaporation process. 



We choose a lattice by dividing space into cells, the size 
of the j th cell being <jj, and we go to the limit <jj — > 0, Vj. 
The wave functional and all results are completely inde- 
pendent of the choice of the Oj . We define the integral in 
the exponent of l|l(J|l as the limit of a Riemann sum over 
the lattice points as the spacing goes to zero. 



dxf'(x)W(r(x),R(x),f(x)) 



CTj— (Jj 



The functional in (|10f) is then a product state over the 
individual shell wave-functions and can be expressed as 



J 



* [r, R, f] = Y[ * j ( Tj , Rj ,/,-) = II Cxp 



-^(fj+i-mrj+UARjJj)} 



(12) 



r 



The energy, Sj, contained between lattice points j and 
j + 1 is related to fj+i — fj by 

fi+i-fi = ?%M +1 -Mj)<? = ^ (13) 



implying that the wave-functional can be expressed in 
the form 



n« 



(14) 



where we have defined Sj = hujj = (Mj + i — Mj)c 2 . Our 
strategy is to solve the Wheeler-DeWitt equation © in- 
dependently at each label, j. The continuum limit would 
be defined by putting together the shell wave-functions 
and taking the limit as the lattice spacing approaches 
zero. We shall see that the wave functions are indepen- 
dent of the original choice of cell size, as it should be. 

Before proceeding further it is necessary to define what 
is meant by a functional derivative when functions are 
defined on a lattice. The defining equations can be un- 
derstood by analogy with the simplest properties of func- 



tional derivatives of the functions J(x 

SJ(y) 
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5J(x) 
dyJ(y) 



S(y-x), 



8J(x 

From these definitions follows 
S 



1 



SJ(x) 



dyJ{y)<t>{y) = <t>( x )- 



(15) 



(16) 



On a lattice we define, for the lattice intervals Xi and Xj, 
SJ(xi) . , \ ,. Si 



5J(xj) 



A{ Xi 



Xj) = lim — 

ffi->0, Vi (Ti 



(17) 



where Xi labels the i lattice site and Sij is just the 
Kronecker 8, equal to zero when the lattice sites Xi and 
Xj are different and one when they are the same, (jj is 
the size of the i th site, which we take to be arbitrarily 
given. Just as 8{y — x) is only defined as an integrand 
in an integral, so A(xj — xj) should also be considered 
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defined only as a summand in a sum over lattice sites. and 
Hence 

5 , . \ SJ(xi) 



(18) 



} j a i J(x i )(p(a 



lim^ .^Vj A (xj - Xj)<p(xi) = <p(xj) 

i 

I 



(19) 



It follows that 



lim 



1 a 



8J(Xj) a-j^O, Vj Uj 8Jj 

where Jj = J{xj). This is compatible with the formal 



(continuum) definition of functional derivative. With this 
rcgularization, the Wheeler-DeWitt equation on the lat- 
(20) tice may now be written for each lattice site, j, as 



lim 



d 2 i 
± 



o\a*c 2 dT? a 2 8R 2 3 



± 



(21) 



Factoring out the term 1 /cr? the equation becomes inde- 
pendent of the choice of cell sizes, 



1 d 2 



± 



d 2 



c 2 dr 2 dR 2 



± 



(fj + i-f 3 ) 2 
4/2 1^1 



^ 3 (t,R,F) = 0. 



(22) 

We see that the limit of vanishing cell size may be taken 
and that the wave-functional is invariant under diffco- 
morphisms in this limit by construction. Upon making 
some obvious replacements, 



c) 2 



Uj ( 1 



* j (R j ,F j )=0, 



(23) 



The upper sign refers to the exterior (R > F) and the 
lower to the interior (R < F). is the time indepen- 
dent part of the wave function. Changing back from the 
coordinate R* to the original coordinate R, the equation 
has the form 



it 



i 

77 



(24) 



where the upper sign refers to the exterior [R > F) and 
the lower sign to the interior (R < F) and It j refers to 
the time independent part of the wave-functional in the 
j th cell. Defining the dimensionless variables Zj = Rj/Fj 



and 7j = FjLOj/c, redefining the time independent part of 
the wave function = yj, and suppressing the subscript 
j throughout, the equation for each shell takes the form 



2 d 2 y z — 1 dy 
dz 2 + 



Yz 2 y = 



(25) 



2 dz 

both in the exterior (z > 1) and in the interior (z < 1). 

IV. SOLUTIONS 

Before examining the solutions of (|25|l we recall some 
results of the WKB approximation. The WKB solutions 
near spatial infinity were shown to be |9j 



exp 



211 



"' l drF'(r±—^) 

c 



(26) 



In addition, finding the WKB approximation functional 
required no regularization. This suggests that, for con- 
sistency, our shell wave-functions should agree with the 
WKB approximation. 

Writing (|25|l in terms of u = l/z and keeping only the 
lowest order term in u we find 



du 2 



2dy 
u du 



7 



y = o, 



(27) 
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which has as a solution 

y ± (z) » ^ ±e ±2^v^+|ln^ 



(28) 



The time dependent wave-function then has the form 



n 



(29) 



As before, the last term in the exponent above is small 
compared with the first two and may be neglected to 
leading order. Then we recover the near-horizon WKB 
solutions of where it was shown that the logarith- 
mic term, which represents the scattering of the wave- 
function in the background geometry, yields Hawking ra- 
diation at the Hawking temperature. In the infinite fu- 
ture, the logarithmic term approaches negative infinity 
in the approach to the horizon, therefore on X + only the 
component with the positive sign is relevant and we have 



We see that the last term is of order h and may be ignored 
when compared with the WKB approximation. Hence 
the wave function we obtain at infinity agrees with the 
results of the WKB approximation at infinity. 

For the asymptotic wave-functions of (|2*§)l to describe 
collapse, we must choose shell wave-functions that cor- 
respond to incoming waves in the infinite past and at 
large distances from the center. These are the wave- 
functions with a positive sign in the exponent and they 
vanish along I + . The time reversed situation, in which 
the wave- functions vanish on X~ but not on X + , are out- 
going modes and represent expanding matter. They are 
described by the wave- functions in l|29|l with the negative 
sign in the exponent. 

Near the horizon, \z — 1| « 0, we define u = z — 1. 
Again, retaining only the lowest orders in u, we find the 
effective equation 



du 2 2u du u 2 
and the two independent solutions 



(30) 



(31) 



hor 



n 



representing outgoing shells, scattered near their hori- 
zons. 

Exact solutions to 125|) may be obtained in a neigh- 
borhood of ordinary points or about non-essentially sin- 
gular points by expanding in a Frobenius series. From 
the point of view of Hawking radiation, the horizon is 
the more interesting of the two non-essential singulari- 
ties of (|25|l . Expanding about the horizon will yield at 
least one solution by Fuch's theorem. Instead of z then, 
it is convenient to use the variable u = z — 1. Assuming 
a solution of the form 



oo 

y (u) = Va„ U fc+ " 



n=0 



we find that the roots 



k± = - (l ± y/1 - 16 7 2 



(34) 



(35) 



of the indicial equation 



where <7j = J 1 — I67 2 . We consider only astrophysical 

black holes, for which M is large. Let 7 j ^> 1, so that 
Oj ~ 4«7j , then the solution may be put in the form 



hor 



^hor 



II' 



-^.( TJ ±4H^-i|+^-H*i-i|) 



(32) 



k\k-- 



7 2 = 



(36) 



differ by a non-integral number and therefore two inde- 
pendent solutions are obtained, each defined by the fol- 
lowing recursion relations, 



[k±(k± -I) + 2 7 2 ]4 



7 2 a„ + [(fc± + n + l)(fc ± + n) + 2 7 2 ]a* 
7 2 + (fc± +n + 2)(k± +n+§) 



a n+2 ~ 7 2 + (fc ±4 " T ' ' 

r 



Each series provides an exact solution for small \z — 1|, combination of the two, 
and can be shown to be convergent for |u| = |z — 1| < 1. 

The general solution at each shell is therefore a linear ~ <^2> , k-++n 

v| './.: "7 2- "./.<•"./ 1 (38) 
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where we have reinserted the shell label, j. The very first 
order in the expansion gives precisely (|33l) and therefore 
Hawking's formula. Because this is the exact quantum 
state of the LTB class of collapse models it also contains 
all the information about the quantum evolution of the 
collapse close to the horizon, up to the scale of validity 
of the canonical theory. We have shown that the first 
order approximation yields the standard Hawking radia- 
tion at the Hawking temperature The higher orders 
represent corrections to Hawking's original formula. 



The measure in the above integral is determined by trans- 
forming to the Schwarzschild Killing time, T, which is 
related to the Tolman-Bondi (proper) time by a stan- 
dard transformation The wave-functionals consid- 
ered must also be expressed in terms of the Schwarzschild 
Killing time. Ignoring the unimportant T-dependent 
part (which will go away on squaring) the Bogoliubov 
coefficient turns out to be 



(43) 



V. HAWKING RADIATION 

We now demonstrate that the leading correction to 
Hawking radiation makes the spectrum non-thermal. 
First we show that the Hawking spectrum itself arises 
from retaining only the n = term and dropping all 
terms of 0(H) in the exponent of Q29J1 . If one further 
assumes that 7> 1, then (|35|) gives k± » 1/4 ± 17, and 
the above shell by shell wave-functions can be written as 
(dropping the shell labels) 



hor 



A hor 



exp 



-iu(r -I In \z - 1|) 

c 



(39) 



and the asymptotic WKB wave- functions in (|29[) are 



« A± exp 



1 ,2F 
- /a' I r ± — v z 



(40) 



Taking the negative sign in the exponent, the function 
ty 00 in (|4*U)) represents a free outgoing dust mode of fre- 
quency u>. We take a complete basis of outgoing modes 
to be given by 



exp 



2F 

c 



(41) 



letting lj' take all possible values between zero and infin- 
ity. 

The states v[/ hor are the scattered modes near the hori- 
zon. We are interested in computing the projection of 
(|3l?)) on the negative frequency modes of the outgoing ba- 
sis. For this purpose, we must consider the inner product 
of the states on a hypcrsurface of constant Schwarzschild 
time. As explained in [jj, the Bogoliubov /3-coefficient of 
interest is given by 



zdz 



1 



This integral can be performed, after making the substi- 
tution s — \fz — 1. Retaining only the lowest order in s 
one gets 



0(u),u)') « 2Fe 



dse^^ s [2s 



sr „l-l-2iJ 



(44) 



After squaring, one obtains the desired transition proba- 
bility as 



\2_o„2 T? 2 kT H 1 



(w,u/)| 2 = 27r 2 F 2 



£ _ 1 



(45) 



which is the Hawking spectrum at the Hawking temper- 
ature, kTn = he 3 /8irGM . Note that M is the total mass 
contained within the radiating shell. 

Apart from neglecting all 0(H) corrections to the wave- 
functionals, the Hawking formula is recovered by employ- 
ing a "near horizon" approximation in which only the 
lowest order in s — y/z — 1 is considered. The latter is 
justified by the fact that only the near horizon region 
contributes significantly to the integral in (|43() through 
the pole at z = 1. 

Consider the correction to this formula, which is ob- 
tained by taking into account the first 0(H) correction. 
Keeping the asymptotic form of the wave-functional the 
same, as given in (|41f) . but including the first 0(K) correc- 
tion to the near horizon wave- functional in (|33|l , we eval- 
uate the Bogoliubov coefficient in l|42() to leading approx- 
imation in s {i.e., near the horizon). With s = \fz — 1, 
this is 



P(w,u') w 2Fe 



dse~ 



\2s] 



(46) 



Although it appears that the r.h.s. of 1)46(1 has the 
same structure as the WKB expression, it is evident that 



(42) \/3(uj,u)')\ 2 is no longer thermal. When 7 > 1 we find 



( W , W ')| 2 «27r 2 F 2 



2c kT H 1 
Flo' e e fcTj7 — 1 
I 



-Re ip 



is 



2nkT H 



(47) 



where ip(z) is the Polygamma function. As the first correction depends on the Hawking temperature (or the 



7 



mass M) and not just the emission frequency, it cannot 
be associated simply with a correction to the density of 
states. Furthermore, it cannot be obtained by modify- 



J 



ing the Hawking temperature. When 7 > 1, the real 
part of the Polygamma function can be approximated by 
Re[V>(— iy)\ ~ my and we find, 



{uu,uu')\ 2 ^2k 2 F 2 



2c kT H 1 



Fu' 



e kT H _ I 
I 



■In 



7rfcT, 



H 



(48) 



which completes the first correction to the WKB approx- 
imation. The correction term is of order 0{hxK). 



VI. 



DISCUSSION 



The correction we have obtained cannot be accounted 
for by modifying the Hawking temperature and we con- 
clude that it renders the radiation non-thermal. The 
breakdown of Hawking's thermal spectrum, when quan- 
tum gravitational effects are accounted for i.e., when we 
go beyond the WKB approximation, is not equivalent to 
a correction to the black hole entropy. Corrections to 
the black hole entropy formula have been computed in 
various approaches before [l^, but they can all be un- 
derstood as relating to a thermal spectrum, simply mod- 
ifying the Hawking temperature. Put in another way, 
we have shown that there are "grey body" factors whose 
origin is quantum gravitational and they are associated 
with the horizon. They are not to be confused with grey 
body factors that are associated with the radiation at 
spatial infinity, which originate in the backscattering of 
the radiation against the classical space-time geometry. 

The thermal character of the Hawking radiation from 
collapsing matter is sometimes invoked to suggest that 
quantum gravitational evolution may not be unitary and 
that information is lost. Our result indicates that this is 
not necessarily the case, at least within the context and 
limitations of the canonical theory. A more complete pic- 



ture will be obtained from a deeper understanding of the 
behavior of solutions to (|25[) . We hope to come back to a 
more developed discussion of the properties of the solu- 
tions to this equation in a future publication. One of the 
issues we expect to discuss will be the nature of the sin- 
gularity at infinity. It is an essential singularity so there 
will be an issue regarding normalizing the wave function. 
The collapsing system we study has been a marginally 
bound system - physically this means that the total vac- 
uum energy is zero. From quantum considerations we 
expect there to be a non-vanishing zero point energy. A 
preliminary estimate of the result of allowing a non-zero 
vacuum energy is that the singularity at infinity will be 
softened but that the general features of the solution dis- 
cussed so far will be unchanged. 

It was especially gratifying to see that the results are 
completely independent of the choice of cell sizes in the 
lattice. 
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